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Abstract 

We present a detailed comparison of a variety of predictions for diffractive 
light vector meson production with the data collected at the HERA collider. 
All our calculations are performed within a dipole model framework and make 
use of different models for the meson light-cone wavefunction. There are no 
free parameters in any of the scenarios we consider. Generally we find good 
agreement with the data using rather simple Gaussian motivated wavefunctions 
in conjunction with dipole cross-sections which have been fitted to other data. 
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1 Introduction 



In a previous paper, Forshaw, Kerley and Shaw (FKS) [Tj reported on a successful at- 
tempt to extract the cross-section for scattering colour dipoles of fixed transverse size 
off protons using both electroproduction and photoproduction 7p total cross-section 
data, together with the constraint provided by the measured ratio of the diffraction 
dissociation cross-section to the total cross-section for real photons. Subsequently, 
the same model has been applied to "diffractive deep inelastic scatterring" (DDIS) 
7*p — > X p [2] and also to deeply virtual Compton scattering (DVCS) j*p — > -yp [3]. 
In both cases, the model was shown to yield predictions in agreement with the data 
[HE] with no adjustable parameters. The model can also be extended to diffractive 
vector meson production: 

J*(q)+p(p)^V(q')+p(p') V = p,4> or J/*, (1) 

where the squared centre of mass energy s = W 2 = (p + qf > <2 2 ,M£. In these 
processes the choice of vector meson, as well as of different photon virtualities, allows 
one to explore contributions from dipoles of different transverse sizes. The process 
also has the advantage that there is a wide range of available data. It ought also to 
provide important information on the poorly known light-cone wavefunctions of the 
the vector mesons. 

The aim of this paper is to confront the predictions of the dipole model with the 
HERA data on p and <j) electroproduction. The predictions for J/^f production are 
best considered in conjunction with an analysis of open charm production, and will be 
discussed elsewhere. We shall focus on the FKS model Q, but we also compare with 
the predictions of two other models: the Golec-Biernat-Wiisthoff (GW) saturation 
model [SJ Ej; and the recent "Colour Glass Condensate" (CGC) model of Iancu, 
Itakura and Munier [Hj. For the meson light-cone wavefunction, we shall consider 
three different ansatze: the Dosch, Gousset, Kulzinger and Pirner (DGKP) [§] model; 
the Nemchik, Nikolaev, Predazzi and Zakharov (NNPZ) ^U] model; and a simple 
"boosted Gaussian" wavefunction, which can be considered as a special case of the 
latter. 

The paper is laid out as follows. In the first two sections we summarise the dipole 
models used and discuss the forms chosen for the vector meson wavefunctions. We 
then compare their predictions with experiment before drawing our conclusions. 



2 The colour dipole model 

In the colour dipole model ^T] , the eigenstates of the scattering (diffraction) operator 
are "colour dipoles", i.e. quark-antiquark pairs of transverse size r in which the 
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quark carries a fraction z of the photon's light-cone momentum*. In the proton's rest 
frame, the formation of the dipole occurs on a timescale far longer than that of its 
interaction with the target proton. Because of this, the forward imaginary amplitude 
for singly diffractive photoprocesses 7p — > Xp is assumed to factorise into a product 
of light-cone wavefunctions associated with the initial and final state particles 7 and 
X and a universal dipole cross-section a(s,r), which contains all the dynamics of the 
interaction of the qq dipole with the target proton. In particular for reaction $1} one 
obtains 




where ^^(r, z) and ^^(r, z) are the light-cone wavefunctions of the photon and 
vector meson respectively. The quark and antiquark helicities are labelled by h and h 
and we have suppressed reference to the meson and photon helicities. The dipole cross- 
section is usually assumed to be flavour independent * and, as implied by our notation, 
"geometric", i.e. for a given s, it is assumed to depend on the transverse dipole size, 
but not the light-cone momentum fraction z. The light-cone wavefunctions do depend 
on the quark flavour, via their charges and masses. Finally, the corresponding real 
part of the amplitude (j2J) is either neglected or, as here, is estimated using analyticity. 

While the photon light-cone wavefunction can be calculated within perturbation the- 
ory, at least for small dipole sizes, the vector meson light-cone wavefunctions are not 
reliably known, and must be obtained from models. This will be discussed in the 
following section. The rest of this section is devoted to the dipole cross-section, for 
which we shall consider three different models*. Since full details are given in the 
original papers, our treatment will be brief. 

2.1 The FKS model 

The FKS model El 0] is a two-component model 

a(s, r) = a so{t (s, r) + (X har d0, r) , (3) 

in which each term has a Regge type energy dependence on the dimensionless energy 
variable r 2 s\ 

^^^(l-^)^ (4) 

<Thard(s, r) = (of r 2 + a% r 6 ) exp(-z/^r) {r 2 s) x » (5) 

* We work in light-cone coordinates x^ = (x + ,x~,x) in the convention where x^ = x°±x 3 . Here 
z = k + / q + where the momentum k of the quark is (fc + , k~ , fc). 

^For the GW model, this is only strictly so at large Q 2 since some flavour dependence enters 
indirectly at small Q 2 through the definitions of x mo d (see below). 

'■For a more general review of phcnomenological dipole models, see for example |12j . 
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This parametric form§ is chosen so that the hard term dominates at small r and goes 
to zero like r 2 as r — > in accordance with ideas of colour transparency, while the 
soft term dominates at larger r ^ 1 fm, with a hadron-like soft pomeron behaviour. 
In addition, to allow for possible confinement effects in the photon wavefunction at 
large r, FKS modified the perturbative wavefunctions ipj^ L (r, z) by multiplying them 
by an adjustable Gaussian enhancement factor: 

\^ L (r,z)\ 2 = \r T , L (r,z)\ 2 f(r) (6) 

where 

l+Bexp(-c 2 (r-R) 2 ) 
J[) l + Bexp{-c 2 R 2 ) " KJ 

This behaviour is qualitatively suggested by an analysis [TH| of the scattering eigen- 
states in a generalised vector dominance model jHj which provides a good description 
of the soft Pomeron contribution to the nucleon structure function F 2 on both protons 
and nuclei The free parameters in both the dipole cross-section and the photon 

wavefunction were then determined by a fit to structure function and real photoab- 
sorption data. The resulting values are given in Table 1. Having been obtained in 
this way, they were then used to predict successfully the cross-sections for other pro- 
cesses which depend solely on the dipole cross-section and the photon wavefunction, 
namely diffractive deep inelastic scattering (DDIS) [2] and Deeply Virtual Compton 
Scattering (DVCS) 0. 



The resulting dipole cross-section is shown Fig. [TJ As can be seen, as s increases the 
dipole cross-section grows most rapidly for small r, where the hard term dominates, 
eventually exceeding the typically hadronic cross-section found for dipoles of large 
r ~ 1 fm. This rise could well be tamed by unitarity or saturation effects |17j . 
However, the authors have argued ^H] that such saturation effects are unlikely to be 
important until the top of the HERA range and beyond, and they are not included 
in the FKS model in its present form. 



2.2 The GW model 

This well-known model jHJ [7j combines the approximate behaviour a — > r 2 f(x) at 
small r together with a phenomenological saturation effect by adopting the attrac- 
tively simple parametric form: 

o — <jq ^1 — exp 

§This form, taken from [3], is actually a simplified form of that used in the QIE], but gives almost 
identical results. 

^For a more recent discussion of the relation between GVD models and the dipole approach, see 



-r 2 Ql 
4(x mod /x ) ; 



(8) 
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Table 1 



As 


0.06 ±0.01 


X H 


0.44 ±0.01 


a% 


30.0 (fixed) 


af 


0.072 ±0.010 


af 


0.027 ±0.007 


af 
vh 


1.89 ±0.03 
3.27 ±0.01 


B 


7.05 ±0.08 


c 2 


0.20 (fixed) 


R 


6.84 ±0.02 
0.08 (fixed) 




1.4 (fixed) 



Table 1: Parameters for the FKS model [3] in appropriate GeV based units. 
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Figure 2: The GW dipole cross-section (left) and CGC dipole cross-section (right) at 
W = 75 GeV for Q 2 = 2 GeV 2 and Q 2 = 20 GeV 2 . The Q 2 -independent FKS dipole 
cross-section at the same energy is shown for comparison. 

Here x mo d is a modified Bjorken variable, 



where rrif is the quark mass and Q = 1 GeV. The 3 free parameters xq, o"o and A were 
successfully fitted to F 2 data. The four-flavour fit which we shall use in this paper 
yielded <To = 29.12 mb, A = 0.277 and xq = 0.41 x 10 -4 . The quark masses are chosen 
to be 0.14 GeV for the light quarks and 1.5 GeV for the charm quark. The model 
is also able to describe F®^ 1 data [7]. A recent refinement of the model takes into 
account corrections due to DGLAP evolution at large Q 2 but these are rather 
small corrections and are not included here. Finally, all these results are obtained 
with a purely perturbative photon wavefunction, which is somewhat enhanced at large 
r- values by the use of a lighter quark mass than that used in the FKS model". 

The resulting behaviour of the dipole cross-section is illustrated and compared to that 
of the FKS model in Figure 121 



The dipole model of Iancu, Itakura and Munier [Hj can be thought of as a development 
of the Golec-Biernat-Wiisthoff saturation model. Though still largely a phenomeno- 
logical parameterisation, the authors do claim that it contains the main features of 
the "color glass condensate" regime, where the gluon densities are high and non-linear 

"See, for example, Fig. 3 of 




(9) 



2.3 The CGC model 
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effects become important. In particular, they take 



rQ s < 2 

rQ s > 2 , (10) 

where the saturation scale Q s = (x /x) x ^ 2 GeV. The coefficients a and b are uniquely 
determined by ensuring continuity of the cross-section and its first derivative at rQ s = 
2. The leading order BFKL equation fixes 7 S = 0.63 and k = 9.9. The coefficient Ao 
is strongly correlated to the definition of the saturation scale and the authors find 
that the quality of fit to F2 data is only weakly dependent upon its value. For a 
fixed value of Ao, there are therefore three parameters which need to be fixed by a 
fit to the data, i.e. Xq, A and R. In this paper, we take N = 0.7 and a light quark 
mass of m q = 140 MeV, for which the fit values are Xq = 2.67 x 10~ 5 , A = 0.253 and 
R = 0.641 fm. 

As for the GW dipole, we compare to the FKS dipole at two values of Q 2 in Figure 

121 



a 



2nR 2 J\fn 



«A ln(l/i) J 



2txR 2 {1 - exp [-a In 2 (6rQ s )]} 



for 
for 



3 Light-cone wavefunctions 



The light-cone wavefunctions ^ h ^(r, z) in the mixed representation (r, z) used in the 
dipole model are obtained from a two dimensional Fourier transform 

*H*{r,z) = J ^e^^M) (11) 

of the momentum space light-cone wavefunctions i S h ^(k, z), where the quark and 
antiquark are in states of definite helicity, h and h respectively. For transversely 
(T) or longitudinally (L) polarised photons, the momentum space light-cone wave 
functions themselves are calculated perturbatively (HIED! (per fermion of charge ee/): 

**>(M) = yg^ (ee/7 ^^§ 4 ,( M ) . (12 ) 

Hereafter, A denotes the polarization state L or T. e x are the polarisation vectors 
and the "scalar" part** of the photon light-cone wavefunction, $ 7 , is given by 

*"< k -' ) = m-w + l + m y (13) 



**This would indeed be the photon light-cone wavefunction in a toy model of scalar quarks and 
photons. 
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For the vector mesons, the simplest approach is to assume the same vector current 
as in the photon case, with an additional (unknown) vertex factor T\(k, z): 



<«(M ^^(7-4)^#>r(M (i4) 



where the scalar part of the meson light-cone wavefunction is given by 

AV 1 -2(1- z)M$ + k 2 +mj v 1 

Different models are defined by specifying these scalar wavefunctions. In practice, it 
is common to choose the same functional form for <3>^ and perhaps allowing the 
numerical parameters to differ. 

It is instructive to consider the longitudinal wavefunctions more explicitly Using the 
polarisation vectors 

and the rules of light-cone perturbation theory given in [20] , it follows that the lon- 
gitudinal photon light-cone wavefunction is 

*Js(M = ^K-m y k , + m) + z{l _ z)Q , - q) ■ (17) 

and that of the vector meson is 

^ (M) _/A^ _{ z(l-z)2M v T(k,z 




mr 



- z(l - z)M% 




On substituting (|T7j) in (JTTJ) the second term of (|T7j) leads to a dipole of vanishing 
size, which does not contribute to the cross-section. This is in accord with gauge 
invariance. The same argument cannot be used to justify the omission of the second 
term in the meson wavefunction ()18|) . since the latter has a k dependence. In practice, 
this term is omitted in the DGKP model 0, but retained in the NNPZ model ^U]- A 
discussion of the gauge invariance issues surrounding this point can be found in 



Before discussing these models more fully, we give the explicit forms for the photon 
wavefunctions in r-space. 

The normalised photon light-cone wavefunctions are jU]: 



and 




Kh^ z ) = \l 4^,-^/2^(1 - z)Q^^ , (19) 



<£ =±) (r,^) = ±y^ee / ^ ± ^(^ T -(l- 2 )^ ± )9 r +m / ^ ± ]^^ , (20) 



where 

e 2 = z(l - z)Q 2 + m) . (21) 

Since the modified Bessel function K (x) decreases exponentially at large x, these 
equations imply that at high Q 2 , the wavefunctions are suppressed for large r unless 
z is close to its end-points values or 1. As can be seen from equations ([Tn|l - (pn]) , 
these end-points are suppressed for the longitudinal but not the transverse case. This 
is the origin of the statement that transverse meson production is more inherently 
non-perturbative than for longitudinal meson production. 

For small r, the perturbative expressions given above are reliable. For large r-values, 
however, confinement corrections are likely to modify the perturbation theory result. 
These larger r-values contribute significantly at low Q 2 , where the wavefunctions are 
sensitive to the non-zero quark masses mf, which prevent the modified Bessel function 
from diverging in the photoproduction limit. For these reasons, the photon light-cone 
wavefunctions at large r are clearly mo del- dependent. 

We now turn back to the meson wavefunctions. These are subject to two constraints. 
The first is the normalisation condition (201 122] 




which embodies the assumption that the meson is composed solely of qq pairs. Note 
that this normalisation is consistent with equation (J2J) and differs by a factor An 
relative to the conventional light-cone normalisation. 

The second constraint comes from the electronic decay width 122] : 

(2n) 2 z(l - z) W- z )Q 2 + k 2 + m 2 )^l- h (k, z)^- h (k, z) , (23) 

h, h 

where the coupling /y of the meson to electromagnetic current can be determined from 
the experimentally measured leptonic width T v ^ e + e - since 3MvT v ^ e + e - = A-Ka 2 m f v . 
We shall prefer to implement the constraint directly in terms of the r-space wave- 
functions. For our purposes, we can write the meson wavefunctions in r-space as 



,V,L /_ , N c „ 1 



n ' n V 4"7T MyZ{Y — Z t 

where V 2 . = -d r + d 2 and 



[z{\ - z)M 2 + 5 x (mj - V r 2 )]0 L (r, z) (24) 



V,T(t=±) 



h.h 



(r,z) = ±1 



^ vg 

47T z(l - Z) 



[ie 1 r (z5 h±rhT -(l-z)5 hTrh ±)d r +m f 5 h±rh ±\(f) T (r J z). 

(25) 
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Note the second term in square brackets which occurs in the longitudinal meson case. 
This is a direct consequence of keeping the second term in (|18|). For the DGKP 
wavefunctions this term is absent, i.e. 5 = 0, whilst NNPZ keep this term, i.e. 5=1. 
In terms of these wavefunctions, (J23J) becomes (assuming that r0y(r, z) — > at r = 
and r = oo) 

f v M v = ^e f f — ^- [z(l-z)M 2 v + 5x(m}-V 2 r )]Mr,z)\ r _ (26) 
71 Jo zy 1 ~ z ) 

and 

fvMv = ~^e f jT 1 RT ^ F [(, 2 + (1 - ,) 2 ) V r 2 - mj] T (r, z)\^ . (27) 



In computing and all other observables involving the meson we in all cases take 
the quark mass to be equal to the light quark mass plus 150 MeV. 



3.1 DGKP meson wavefunction 



In the DGKP approach [§], the r and z dependence of the wavefunction are assumed 
to factorise*. Specifically, the scalar wavefunction is given byt 



<f> v x (r,z) = G(r)f x (z)z(l-z) (28) 
A Gaussian dependence on r is assumed, that is 

G{r) = TT^-e"^ , (29) 
N c t /My 

and f\(z) is given by the Bauer-Stech-Wirbel model |24j : 



h{z) = N x ^z(l - z)e " . (30) 



Setting 5 = (recall that this is equivalent to neglecting the second term in (j!8j) ) in 
(PH) results in 

*J£(r, z) = ,(1 - *) ^-^/^-i^ , (31) 

*Note that the theoretical analysis of Halperin and Zhitnitsky shows that such a factorising 
ansatz must break down at the end-points of z. However, since the latter are suppressed in the 
DGKP wavefunction, this has no practical consequence. 

^DGKP do not actually include the factor z(l — z) in the scalar wavefunction. This is because 
they define the scalar wavefunction to be the r.h.s of l|15|) divided by z(\ — z). 
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for the DGKP longitudinal meson light-cone wavefunction. if is the effective charge 
arising from the sum over quark flavours in the meson: ej = l/V2,l/3 and 2/3 for 
the p, <p and J/^ respectively. Similarly, the DGKP transverse meson light-cone 
wavefunctions can be written as 



± 



lustre 



±ie r 



'I 



X 



Vnfv 



+ ^5 

My 

h{z)e 



h±,h± 



(32) 



The normalisation condition (1221) on the DGKP wavefunction leads to the relations, 



with 



and 



I L = I dzz\\ - zf f L (z) 
'o 



1 [z 2 + {l-zfU + m) 

dz — "/tW 



o Ml 



(33) 

(34) 
(35) 



The leptonic decay width constraints (|26j) and (|27|) on the DGKP wavefunction yield 

- 1 2[z 2 + (1 - z) 2 ]uj 2 + m 2 
dz- 



dz z(l - z)f L {z) 



2M v z(l - z) 



-h(z) ■ 



(36) 



'o Jo 
The values of u>\ and A/a are found by solving (|3^j) and (|3Tj) simultaneously, and are 
given in Tabled The values of the decay constants used are the central experimental 
values |2EI: /p = 0.153 ± 0.004 GeV, and fa = 0.079 ± 0.001 GeV. 

Table 2 



DGKP 


ml 




M L 


N't 


P 


0.331 


0.206,0.218 


15.091 


5.573,8.682 





0.368 


0.244,0.262 


15.703 


5.689,8.000 



Table 2: Parameters and normalisations of the DGKP light-cone wavefunctions in 
appropriate GeV based units. For the transverse case, the first and second values are 
for the FKS and GW quark masses respectively. 



The resulting behaviour of the p wavefunctions is shown for the case of the FKS quark 
masses in Figures 01 As is clear from equations ()28j) - (j32j) . the wavefunctions peak 
at z = 0.5 and r = 0, and go rapidly to zero as z — > 0, 1 and r — > oo, so that large 
dipoles are suppressed. From the figures, we see that the transverse wavefunction has 
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0.02 



0.08 



0.06 



0.04 








Figure 3: The p- wavefunctions |\1/ L | 2 (left) and |\I/ T | 2 (right) in the DGKP model 
with the quark mass used in the FKS dipole model. Note the different scales for the 
ordinate. 

a broader distribution than the longitudinal wavefunction. The <fi wavefunctions are 
qualitatively similar to, but slightly more sharply peaked than, the p wavefunctions. 

The GW model uses a much smaller value for the light quark masses than the FKS 
model, as we saw in Section 2.3. We might expect this to have a striking effect on 
the transverse wavefunction of the p, since the transverse wavefunction ()32)1 vanishes 
at r = for zero quark masses, while the longitudinal wavefunction ()31jl does not. 
The transverse DGKP wavefunction with the light quark mass used in the GW dipole 
model is shown in Figure 0] for the p. As can be seen, the smaller quark mass decreases 
the wavefunction at the origin and shifts the peak to slightly larger r. 

3.2 Boosted wavefunctions 

In this approach, the scalar part of the wavefunction is obtained by taking a given 
wavefunction in the meson rest frame. This is then "boosted" into a light-cone wave- 
function using the Brodsky-Huang-Lepage prescription, in which the expressions for 
the off-shellness in the centre-of-mass and light-cone frames are equated (or equiv- 
alently, the expressions for the invariant mass of the qq pair in the centre-of-mass and 
light-cone frames are equated |27j). 

The simplest version of this approach assumes a simple Gaussian wavefunction in the 
meson rest frame. Alternatively NNPZ ^U] have supplemented this by adding a hard 
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Figure 4: The p-wavefunction |\1/ T | 2 in the DGKP model with the quark mass used 
in the GW dipole model. 



"Coulomb" contribution in the hope of improving the description of the rest frame 
wavefunction at small r. We refer to ^0] for details of this procedure. Here we simply 
state the result, which is that the NNPZ meson light-cone wavefunctions are given 
by equations (j2H 123 with 5 = 1, where the scalar wavefunctions (f>\ (r, z) are taken 
to be a sum of a soft (Gaussian in the rest frame) part and a hard (Coulomb) part: 

, / / / m2 f R2 ^ ( 2z(l-z)r 2 \ 

Ur,z) = A/a Uz(1 - z)V2^exp ( - j exp 1 R2 ] j 

(m 2 R 2 \ . a 3 (r) \ 

x ex P (^-J + 16 ^ A(r^)g3(r, g ) r ^ l(i4(r ' " )r/5(r ' " )} J • (37) 

Here 



A(r,z) = \l + , \ f - 4C 2 a 2 (r)m 2 , (38) 
y z(l — 2) J 

5(r,z) = -£±L, (39) 

a(r) = (40) 

4m/a s (r) 

is a running Bohr radius. The strong coupling a s is chosen to run according to the 
prescription 



and 



4tt 

a s = a for r > r s and a s (r) = — for r < r s (41) 

A, log (l/(A2 CD r 2 )) 

where r s = 0.42 fm, a = 0.8, A QCD = 0.2 GeV and (3 = (33 - 2n/)/3. Apart from 
the normalisation constants A/"a , these wavefunctions depend on two free parameters 
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which are independent of the meson helicity: a "radius" parameter R and a parameter 
C, introduced to control the transition between the hard Coulomb-like interaction and 
the soft confining interaction in the rest-frame wavefunction. The case of a simple 
Gaussian wavefunction in the rest-frame, which we will refer to as a "boosted Gaussian 
wavefunction" can be obtained by simply setting C = (we still choose 5 = 1 when 
considering this wavefunction). 

At this point we comment on two issues associated with the behaviour of the Coulomb 
part of the scalar wavefunction for small r. The first is that (J37|) diverges logarith- 
mically as r — > at z = 0.5. This divergence is however regulated in observables. 
For example, while the resulting squared wavefunctions l^/^^)) 2 exhibit a narrow, 
singular peak at r = 0, z = 0.5, the quantity r|\& y ( A )| 2 which enters the normalisation 
condition (|22|) is zero in this limit. Nevertheless, this singular behaviour does have 
a (finite) effect when computing the meson decay constant which depends upon the 
behaviour of the wavefunction at r = 0. The second issue is that when the scalar 
wavefunction (|37j) is substituted into equations - (J27J), the derivatives in r give 
rise to inverse power divergences at r = when acting upon the running coupling 
()41jl . However, these divergences occur solely in terms which are strictly higher order 
in a s . Henceforth we discard these higher order terms, which is equivalent to omitting 
all derivatives of a s (r) with respect to r when differentiating the scalar wavefunction 
(137)1 . We stress that none of these issues arise when using the boosted Gaussian 
wavefunction. 

It remains to determine the various constants. NNPZ determined both R and C 
by using a standard variational procedure for the initial centre-of-mass wavefunction 
using a non-relativistic potential. They then checked that the resulting predictions 
(12311 were in reasonable accord with the observed leptonic decay widths*. Here we 
follow a slightly modified procedure, since we want to be able to easily adjust the 
quark masses to those assumed in the various dipole models. Specifically, we fixed C 
at the value chosen by NNPZ and vary the value of R to give approximate agreement 
with the decay width constraints (J23J). In practice, we found it adequate to use the 
same i?-value for both the FKS and GW mass choices. The resulting values of R 
and C, with the associated values of the normalisation constants, are shown in Table 
H3 and are not very different from the original parameters of NNPZ. In addition we 
show the results for the boosted Gaussian case (C = 0) in Table EJ 

The behaviour of the resulting p wavefunctions is shown in Figures El with the FKS 
choice of quark masses. The divergence at r = 0, z = 0.5 is not visible since we do 
not plot down to r = 0. Like the DGKP wavefunctions, they peak at z = 0.5 and 
r = 0, and go rapidly to zero as z — ► 0, 1 and r — * oo. However, on comparing these 

*Note that a shortcoming of this model is that equations (|23J) - H25(1 give slightly different pre- 
dictions for the decay constant for the case of transverse and longitudinal meson helicities, because 
NNPZ use helicity independent values of R, C. 
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Table 3 



NNPZ 


R 2 


C 


M L 


M T 


ML) 


fv{T) 


P 


25.0 


0.36 


0.0123,0.0121 


0.0125,0.0137 


0.143,0.147 


0.157,0.109 


<\> 


18.0 


0.53 


0.0122,0.0124 


0.0124,0.0131 


0.078,0.078 


0.087,0.067 



Table 3: Parameters of the NNPZ light-cone wavefunctions in appropriate GeV based 
units, together with the resulting values for the decay constants fy. Where two 
values are given, the first and second values are for the FKS and GW quark masses 
respectively. 



Table 4 



Gaussian 


R 2 


M L 


N't 


ML) 


MT) 


P 


12.3 


0.0213,0.0244 


0.0221,0.0259 


0.153,0.161 


0.203,0.192 





10.0 


0.0214,0.0243 


0.0219,0.0251 


0.075,0.079 


0.095,0.088 



Table 4: Parameters of the boosted Gaussian wavefunction in appropriate GeV based 
units, together with the resulting values for the decay constants fy. Where two 
values are given, the first and second values are for the FKS and GW quark masses 
respectively. 

figures with each other, and with Figures El two clear differences emerge. 

Firstly, the peak in z is less sharp in the boosted Gaussian case than in the DGKP 
and NNPZ cases. 

Secondly, the large difference between the longitudinal and transverse wavefunctions 
found in the DGKP case is much less marked in the NNPZ and boosted Gaussian 
wavefunctions. In both cases, the peak in the transverse wavefunction is still broader 
than that in the longitudinal wavefuction, but it is a small effect in comparison with 
the DGKP case. This presumably reflects the fact that in the NNPZ and boosted 
Gaussian wavefunctions, the parameter R has the same value for both helicities, since 
both wavefunctions are generated from the same non-relativistic wavefunction by the 
Brodsky-Huang-Lepage procedure. 

Figures El show the wavefunctions in the case of the boosted Gaussian wavefunction. 
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Figure 5: The p-wavef unctions \ty L \ 2 (left) and |\I/ T | 2 (right) in the NNPZ model with 
the quark mass used in the FKS dipole model. 
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4 Real parts, slope parameters and cross-sections 



We now have all the ingredients required to calculate the absorptive parts (0) of the 
forward amplitudes for vector meson production. For the case of the NNPZ wave- 
functions, we substitute the photon wavefunctions^ ()19j). (J2"U|) and the vector meson 
wavefunctions (f24*|) . (|23|) with 5 = 1 into ((2)). Summing over the quark/antiquark 
helicities and averaging over the transverse polarisation states of the photon, we ob- 
tain^: 

W^pz = ^{^^ ^ jd 2 ra{s,r) J q dz 

{{m) + z(l - z)M$]K (er)<l> L (r, z) - K (er) V r 2 )0 L (r, z)} (42) 

T _ N c e f V4ira cm f 2 f 1 dz 



(2tt) 2 J K,J J z(l-z) 

[(z 2 + (1 - z) 2 )d r K (er)d r Mr, z) + mJ^ o (er)0 T (r, z)] , (43) 

where 4>L,T(r,z) are given by (pT7|). Similarly, using the DGKP wavefunctions (}3~T]) 
and (jS2J), we obtain |Hj 

W^ GKP = 8 J d 2 ra(s,r) j\ z J°^f v z(l - z)] L {z)e~^ l 2 2z{\ - z)QK (er) , 

(44) 



Qfm^4Q GK p = s I d z ra{s,r) / dz 



(s,r) j\z^f v f T (z)e-^ 2 / 2 



^T er u2 , n ^2i ^ , m / 



^[^ + (l-^]ir 1 (er) + ^X (er)) . 



So far, we have focussed on the imaginary amplitude. Taking into account the real 
part contribution, the differential cross-section is given by 



da 1 ' 



T.L 



dt 



1 



t=o 16%s ' 



5nU J 'Y(l + /n (46) 



where /3 is the ratio of real to imaginary parts of the amplitude. It is most straight- 
forward to reconstruct the real part of the amplitude in the FKS dipole model, where 
the dipole cross-section (JHJ), and hence the amplitude (J2J), is given as the sum of hard 
and soft Regge pole terms. In this case, the real part is given by 



$teA FKS = -SmAoft cot {^^- J - 9taoA»d cot , (47) 

§For the FKS case, we must also include the effect of the enhancemect factor (jTJ). 
^ These results reduce to equation (8) of ^U] if, in l|42|) . we integrate the second term by parts 
assuming an r-independent dipole cross-section. 
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where ols,h = 1 + ^s,h = 1-06, 1.44 and Siiu4. soft and Snv4 hard are the contributions 
from the soft and hard Pomeron pieces of FKS dipole cross-section (J3J) respectively. 

It follows from (|47jl that, in the FKS model, (3 lies between 0.09 and 0.83, corre- 
sponding to pure soft and hard pomeron dominance respectively; within this range, 
the value of (3 reflects directly the relative importance of the hard pomeron. Other 
things being equal, (3 will therefore increase with increasing energy, because the hard 
term in the dipole cross-section increases more rapidly with energy than the soft term. 
It will also increase with Q 2 , because the hard term is dominant for small dipoles, 
which are increasingly explored as Q 2 increases. These features are illustrated in 
Figure which shows the values of (3 obtained as functions of W and Q 2 obtained 
using the DGKP wavefunctions. Here we also see that (3 is larger for the longitudi- 
nal than for the transverse case, reflecting sharper peaking of the longitudinal vector 
meson wavefunctions. Similar results for the real parts are obtained using the NNPZ 
wavefunctions. 




W[GeV] Q 2 [GeV 2 ] 



Figure 7: Left: ^-dependence of the ratio (3 of the real to imaginary part of the 
forward amplitude in the FKS model, using the DGKP wavefunctions, for Q 2 = 10 
GeV 2 . Right: Q 2 -dependence of (3 in the FKS model, using the DGKP wavefunctions, 
for W = 75 GeV. 

From the above figures, we see that while the corrections from the real parts in the 
cross-section formulae (j4*U|l are clearly significant in some kinematic ranges, they are 
nowhere dominant. Because of this and because the ratio (3 is expected to be similar 
in the different models" , we shall use the estimates (J47)) of the ratio (3 obtained in the 
FKS model in both dipole models. 

"This was confirmed explicitly for two distinct dipole models in PJ. 
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Assuming the usual exponential ansatz for the t dependence, the total cross-sections 
are given by 

da T > L /dt\ t=0 



(? T,L {l*P -> Vp) 



B 



(4£ 



Unfortunately, the values of the slope parameter B are not very accurately measured. 
Here we use a parametrisation 



B = 0.60 



11 



(Q 2 + M\ 



+ 1 



v) 



(49) 



obtained from a fit to experimental data by Mellado j29j, and used in their analysis 
of the predictions of the GW dipole model by Caldwell and Soares jSH] • The resulting 
values for B are shown in Figure El where we also show the results of an alternative 
parameterisation of Kreisel jH5 to illustrate the range of values that can be obtained 
from different fits to the data. When comparing the predictions of (}4*%j) for the vector 
meson cross-sections with data, it is important to bear in mind that this uncertainty 
in the input value of the slope parameter can easily introduce errors up to of order 
30% or so and that, within this range, this error may be Q 2 dependent. 




20 30 

Q 2 [GeV 2 ] 



Figure 8: Left: Parametrisation of the slope B by Mellado j29] for the p and 
Parametrisation of the slope B by Kreisel [3^ f° r the p and </>. 



Right: 



Finally, on comparing with experimental data, we show 

a TOT = a T + ea L with e = q gg 
in all our plots, although the HERA data range from 0.96 to 1.00. 



5 Results 



In this section we will compare the predictions of the three dipole models with the 
experimental data for our three choices of vector meson wavefunction, without any 
adjustment of parameters. Before doing so, however, we emphasize again that the 
uncertainty in the slope parameter can give rise to Q 2 -dependent errors of up to 30% 
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in the cross-section, which should be borne in mind when comparing with experiment. 
This effect should, hopefully, be much smaller in the ratio R. However, it ought not 
to be completely negligible since the longitudinal and transverse cross-sections will 
presumably have slightly different slope parameters, as to some degree they explore 
dipoles of different sizes. 

5.1 p meson production 

The predictions of the FKS model for the Q 2 dependence of the cross-section and for 
the longitudinal to transverse ratio R are shown in Figure |H1 for the three different 
wavefunctions considered. The equivalent plots for the GW and CGC model are in 
Figures EH] and [TT] respectively. 

The ZEUS total cross-section data are at the following centre-of-mass energies: 0.27 
GeV 2 < Q 2 < 0.69 GeV 2 , W = 51.1 GeV, otherwise W = 75 GeV. The HI data on 
the longitudinal to transverse ratio are, for Q 2 = 9.8, 18.25 GeV 2 , in the range 40 
GeV < W < 140 GeV, otherwise they are at W = 75 GeV. Similarly the ZEUS data 
are, for Q 2 = 0.33 GeV 2 , at W = 67 GeV or for Q 2 = 0.62 GeV 2 at W = 47 GeV. 
For the theory curves we always take W = 75 GeV. 

For all three dipole models, the data favour the boosted Gaussian wavefunction, whilst 
the DGKP wavefunction produces reasonable agreement for FKS and is rather less 
satisfactory for the GW and CGC models. The NNPZ wavefunction is well below 
the total cross-section data for all three models and, for the GW and CGC models, 
in disagreement with the data on the longitudinal to transverse ratio. However, 
as noted in section 3.2, although the spurious singularity in this wavefunction does 
not contribute directly to the predicted cross-sections, it influences them indirectly 
because it influences the value of the radial parameter R deduced from the decay 
width. In what follows, we will therefore focus on the DGKP and boosted Gaussian 
wavefunctions. 

Figures ^1 EH an d El show the W dependence at fixed values of Q 2 for the p me- 
son. Apart from the normalisation, the jy-dependence is good in all cases. The 
normalisation is best described by the boosted Gaussian wavefunction. With this 
wavefunction, the FKS model is in reasonably good agreement with the data, except 
for the ZEUS data at very low Q 2 ; while the GW and CGC models give reasonable 
agreement everywhere. 

This last comment contrasts somewhat with the work of Caldwell and Soares |3()j . 
who have already presented predictions for the GW model using a boosted Gaussian 
wavefunction. These authors also found good agreement with the data on the longi- 
tudinal to transverse ratio and with the W dependence at fixed Q 2 apart from the 
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normalisation. However their results for the normalisation, and hence, implicitly, for 
the Q 2 -dependence of the production cross-section, were relatively poor. However, 
these authors did not implement the leptonic decay width constraint but fixed the 
radial parameter R in the boosted Gaussian by requiring that the exponential in R 
of the wavefunction gives a value of 1/e when the qq invariant mass is equal to the 
meson's mass. This yielded R 2 = 15.5 GeV" 2 for p and R 2 = 8.3 GeV" 2 for <p com- 
pared to our values shown in Table 4. In addition, they neglected real parts, which 
can result in a 20% reduction in the cross-section for large Q 2 . 




Q 2 [GeV 2 ] Q 2 [GeV 2 ] 

Figure 9: The Q 2 dependence of left the total cross-section and right the longitudinal 
to transverse cross-section ratio for p production at W = 75 GeV in the FKS model 
using the three different meson wavef unctions. Data from left [32] and right |3*2*l I33j . 

5.2 (f) meson production 

The corresponding plots to those in the previous section are repeated for the <f> meson 
in Figures El El an d El The HI total cross-section data are at the following centre- 
of-mass energies: for Q 2 = 7.5,8.3,12.5,14.6,17.3 GeV 2 , W = 100 GeV, otherwise 
W = 75 GeV. Similarly for the ZEUS data, for Q 2 = 0.0,8.2,14.7 GeV 2 , W = 
70, 94, 99 GeV respectively. The HI data on the longitudinal to transverse ratio are 
in the range 40 GeV < W < 130 GeV. Similarly the ZEUS data are in the range 25 
GeV < W < 120 GeV. For the theory curves we always take W = 90 GeV. 

The situation is rather similar, as one might expect, to that of the p, i.e. all three 
dipole models tend to do rather well with either of the DGKP or boosted Gaussian 
wavef unctions, while the NNPZ wavefunction is less satisfactory. 
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Figure 10: The Q 2 dependence of left the total cross-section and right the longitudinal 
to transverse cross-section ratio for p production at W = 75 GeV in the GW model 
using the three different meson wavefunctions. Data from left jH2j and right |82| 133]. 




Figure 11: The Q 2 dependence of left the total cross-section and right the longitudinal 
to transverse cross-section ratio for p production at W = 75 GeV in the CGC model 
using the three different meson wavefunctions. Data from left |H21 and right jH2HHBj- 
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W [GeV] W [GeV] 

Figure 12: The W dependence of the total cross-section for p production at various 
values of Q 2 . We use the FKS dipole model and compare the boosted Gaussian and 
DGKP wavefunctions. Data from left jHlj and right |32j . 




W[GeV] W [GeV] 



Figure 13: The W dependence of the total cross-section for p production at various 
values of Q 2 . We use the GW dipole model and compare the boosted Gaussian and 
DGKP wavefunctions. Data from left and right |32j . 
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Figure 14: The W dependence of the total cross-section for p production at various 
values of Q 2 . We use the CGC dipole model and compare the boosted Gaussian and 
DGKP wavefunctions. Data from left and right |H2j . 




Figure 15: The Q 2 dependence of left the total cross-section and right the longitudinal 
to transverse cross-section ratio for production at W = 90 GeV in the FKS model 
using the three different meson wavefunctions. Data from left jS^l EH] and right 

|23 EH. 



24 




Q 2 [GeV 2 ' 



Figure 16: The Q 2 dependence of left the total cross-section and right the longitudinal 
to transverse cross-section ratio for production at W = 90 GeV in the GW model 
using the three different wavefunctions. Data from left jHSl IHEj and right jH3 1 



o 10 




Q 2 [GeV 2 ' 



Figure 17: The Q 2 dependence of left the total cross-section and right the longitudinal 
to transverse cross-section ratio for <fi production at W = 90 GeV in the CGC model 
using the three different wavefunctions. Data from left jHSHSEj and right pTT^ IH8j . 
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6 Conclusions 



We have performed a detailed study, comparing the predictions on p and meson 
electroproduction arising from three different models of the meson wavefunction in 
combination with three different models for the fundamental dipole cross-section. Our 
results are broadly encouraging and support the use of the dipole approach. 

The data can be explained rather well using the dipole model of Forshaw, Kerley and 
Shaw, or those of Golec-Biernat and Wiisthoff and of Iancu, Itakura and Munier, in 
conjunction with either the boosted Gaussian or DGKP meson wavefunction. Cer- 
tainly, we anticipate that excellent agreement could be obtained if one decided to tune 
the meson wavefunctions. Note that agreement extends to the ratio of longitudinal 
to transverse meson production. The NNPZ wavefunction, which has, as noted, an 
unphysical singularity at z = 0.5, r = 0, is not so successful. 

For the future, it is clear that one could use the high quality data from HERA to 
constrain the meson wavefunctions provided the dipole cross-section is sufficiently 
constrained. 
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